Cytoskeleton confinement of red blood cell membrane fluctuations 
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We analyze both the static and dynamic fluctuation spectrum of the red-blood cell in a unified 
manner, using a simple model of the composite membrane. In this model, the two-dimensional 
spectrin network that forms the cytoskeleton is treated as a rigid shell which is located at some 
constant average separation from the lipid bilayer. The cytoskeleton thereby confines both the 
static and dynamic fluctuations of the lipid bilayer. The predictions of the model account for the 
wavevector and frequency dependence of the experimental data. The observed amplitude of the 
thermal fluctuations is related to effects of ATP-driven fluctuations. 



(N 
(N 



o 

CO 
C3 



T3 

C 

O 

o 



> 

y—i 

in 
o 

(N 
O 
-i— » 



i 

-a 
c 

o 
o 



X 



PACS numbers: 

A long-standing problem in the study of red blood 
cell (RBC) structure is the simultaneous softness of its 
membrane observed by thermal fluctuations Q, and the 
strong shear elasticity found in static deformation ex- 
periments, such as micropipette aspiration Q and elec- 
trodeformation S . The membrane itself is a composite 
structure Q with an outer, gel-like extracellular network 
of long sugar molecules (thought to be irrelevant to the 
structural strength), a mixed lipid/protein bilayer and 
an attached, intracellular network. Previous theoretical 
models of this membrane treated it as a single, polymer- 
ized network with the combined curvature bending mod- 
ulus of the lipid bilayer k and the shear rigidity of the cy- 
toskeleton fi S, Such models were successful in describ- 
ing the response of the membrane in static deformation 
experiments, which give [i ~ 10~ 5 — 10~ 6 J/m 2 ^ |^, |J. 
However, comparing these models to the fluctuation data, 
leads to the conclusion that the membrane behaves as if 
the shear modulus vanishes ^ ~ f^, [?]] ■ This surprising 
conclusion comes from the shape and amplitude analysis 
of the longest wavelength shape fluctuations [Q, || . Var- 
ious ideas have been raised in order to account for this 
observation, the main suggestion being that ATP-driven 
structural rearrangement in the spectrin network |l|. || 
relaxes the shear- like deformations. 

The previous studies were concerned with the shape 
fluctuations of largest wavelength. Here, we focus on the 
fluctuation spectrum, at the length-scales of 1 — 0.01/im, 
where the effect of the cytoskeleton is observed ||. The 
important question is to what extent the cytoskeleton ef- 
fects are distinguishable from the fluctuations of a free, 
closed bilayer. We show that for a consistent description 
of both the static and the dynamic fluctuation spectrum 
we must include the confining effects of the cytoskele- 
ton. A simple model where the only significant effect 
of the spectrin network is to confine the lipid bilayer 
membrane, consistently describes both the spatial and 
temporal spectra of the thermal fluctuations of the RBC 
membrane. We discuss how active processes that effec- 
tively increase the temperature of the membrane, can be 
accounted for within this description. 

The curvature bending modulus k of the lipid bilayer 



]io| is deduced from measurements of the amplitude 
of thermal fluctuations at the smallest measured wave- 
lengths | k ~ 2ilx 1(T 20 J. The RBC cytoskele- 
ton is a two-dimensional, roughly triangular, network 
of spectrin proteins |n]], that is attached to the lipid 
bilayer at the nodes and at additional, random sites 
along the spectrin polymers. The cytoskeleton is well 
described as a network of entropic springs, of length 
I ~ 80nm @ [[1, with an effective spring constant: 
/j ~ 4 x 10~ 6 J/m 2 , which is close to the measured 
static shear modulus: ~ 6 x 10~ 6 J/m 2 0]. Compared 
with the lipid bilayer bending modulus k, the curvature 
bending modulus of the cytoskeleton is ncglig ible @, fL3): 
K-cyto ~ l^w 2 ~ 10~ 21 — 10~ 22 J for a cytoskeleton thick- 
ness of w ~ 100 - 500A. 

We now begin by analyzing the measured spatial cor- 
relations of the fluctuations || , and describe the effects 
of the cytoskeleton on the bilayer in terms of continuum 
mechanics. This is feasible since the cytoskeleton forms 
a rather open mesh that is attached to the bilayer at dis- 
crete points with a contact area that is small (~ lnm) 
compared to the inter-node distance (~ lOOnm). In a 
coarse grained picture, we can describe the thermal fluc- 
tuations of the bilayer using a continuum model of the 
free energy functional 



F 



dS 



h h2 



(1) 



that includes the bending and effective surface tension 
energy of the bilayer (see for example |l4j ) in terms of the 
normal displacement h. The surface tension coefficient 
a arises from the constraint of constant surface area of 
the bilayer in the closed geometry of the RBC |l0|, fl5| : 
(j = an/R 2 ~ lx 10 _9 J/m 2 (taking R ~ 4/Ltm as the 
RBC radius, a = 1). The constant can vary awl — 100. 

The last term in (^) is mathematically equivalent to 
a Lagrange multiplier that constrains the mean square 
amplitude of bilayer fluctuations to be equal to: d 2 = 
/cbT/8(7k) 1 / 2 , for an infinite bilayer. This term describes 
the effect of the cytoskeleton on the bilayer through a 
harmonic potential that maintains an average separa- 
tion d (of order w) between the lipid bilayer and the 
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cytoskclcton |TJ|, here treated as a separated, infinitely 
rigid shell, that does not participate in the thermal fluc- 
tuations. The discrete contacts that maintain the con- 
stant average separation are not specifically described in 
this continuum model; in a coarse-grained picture, these 
contacts are the physical origin of the constraint (po- 
tential) that determines the average membrane-spectrin 
network separation. 

The attachment of the cytoskeleton to the bilayer also 
causes stretching and undulations of the bilayer |L6| , 
partly due to steric repulsion between the spectrin and 
the bilayer around the point of attachment |17j . Balanc- 
ing the spectrin stretching with a local curvature of the 
bilayer, results in a membrane with undulations of wave- 
length @, p| L ~ yfnjji ~ 100 - 200nm and ampli- 
tude ~ 10nm(~ w) |l], ^3). In our confining-shell model 
this length-scale is related to the potential-induced per- 
sistence length of the bilayer £0 = (^/t) 1 ^ 4 ~ L (m), i.e. 
the wavelength below which the bilayer is freely fluctu- 
ating. 

We now calculate the spatial correlations for a two- 
dimensional, flat bilayer, since for all but the largest 
wavelengths A, the surface of the RBC is relatively flat: 
50nm < A < lfj/m < R ~ 4/im. We determine the values 
of a and 7 by fitting to the experimental data. From 
Eq. ([!]) the equal-time (static) correlations of the normal 
deflections of the bilayer can be written j|, |l4| 
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In the inset of Fig.l we plot the measured value of 
Kg Q in the form (k q /k — 1) _1 as a function of the 
normalized wavevector (qd) 4 (where d is determined by 
fitting the data to obtain 7 oc d~ A ). From the lin- 
ear slope in the limit of q — > we find the values of 
the parameter 7 = 7.5,1.0 x 10 7 J/m 4 , for the two cells 
measured. These values correspond to mean amplitudes 
d ~ 200,3501 and £ = («/7) 1/4 = 130,220nm respec- 
tively. At larger values of q there is a noticeable de- 
viation from a straight line, which arises from the ef- 
fective surface tension a ~ 7.7,2.8 x 10~ 7 J/m 2 for the 
two cells. Note that surface tension alone, without the 
confining effect of the cytoskeleton (i.e. 7 = 0), does 
not fit the data (dash-dot line, Fig.l inset). In addition, 
there is a rather abrupt change at the crossover wavevec- 
tor qo = l/£o (indicated by the vertical dashed lines in 
Fig.l), above which the data are better described using 
a ~ (To ~ 1.4 x 10~ 9 J/m 2 (solid lines in Fig.l). 

The measured surface tension is consistent with: a ~ 
K /£o ( 0I or der ~ ///lO). This expression gives the ef- 
fect of bilayer shape constraint due to the static undu- 
lations of lateral size £0 x £0, described above. Indeed 
at length-scales shorter then £0, there is no stretching of 
the cytoskeleton (Fig.l). Note that the effective surface 
tension of a closed bilayer is a very sensitive function of 
the excess area of the bilayer [JloL [l|], which is affected 
by the induced undulations. The spread in the measured 



parameters may be due to natural variations in the cy- 
toskeleton network of normal RBC cells. 

There is a qualitative difference in power-law of the 
wavevector dependence of n q for RBC and empty vesi- 
cles p5J[ . The vesicles are well described (Fig. 2) by equa- 
tions (0,||) with 7 = 0, and an effective surface tension: 
Uvesicie ~ k/R 2 ~ 2x 10 _10 J/m 2 , where R ~ 27/xm 
and k ~ 1.3 x 10 _19 J are the vesicle radius and bending 
modulus respectively. Both the RBC's and vesicles data 
collapse when the wavelength is scaled by the r.m.s am- 
plitude d (Fig. 2). For the vesicles of diameter ~ 50/im 
Jl5| the r.m.s. thermal amplitudes are d ~ 1 — 1.5/xm 
(note that here d is not related to confinement). The 
good scaling of the data indicates that there is indeed 
a single important length-scale in the problem, namely 
the persistence length £0 , that determines all the param- 
eters appearing in the free energy (7 and a) and the r.m.s 
amplitude d. 

We now use the same simple model of spectrin confine- 
ment of the bilayer to describe the temporal correlations 
of the membrane fluctuations. The shape fluctuations 
of the RBC membrane are driven by both thermal and 
metabolic energies. The active fluctuations have a fre- 
quency spectrum that is confined to the range 0.3-lHz 
pOj . For higher frequencies, our analysis shows that the 
active processes can be accounted for by an increase in 
the effective temperature of the fluctuations pO, 21]. The 
temporal height-height correlation function [j], 22 1 for a 
flat bilayer at a distance D u from a rigid wall, is 
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where n q is given in (^). The hydrodynamic interac- 
tion (Oseen interaction kernel p3| ) has a modified form: 
A(q) = (1/4 W ) [1 - (1 + 2qD LJ )e- 2 ^] (A f (q) = l/4 V q 
for a free bilayer), so that the relaxation frequency uj(q) 
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where 77 ~ 3i] wa t e r is some average viscosity of the cyto- 
plasm and external solution. In the limit of short wave- 
lengths (q —>■ 00) we recover the free bilayer frequency: 
uj(q) — > Kq 3 /Arj. 

The mean square amplitude of the normal fluctuations, 
as a function of frequency u, is the Fourier transform 
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For a free bilayer this expression (|^) gives an anoma- 



lous frequency dependence fl22j: d(u) 



5/3. We 



integrate the expression (p|) numerically in the range 
-k/R < q < n/a (a ~ 50A), and compare with the ex- 
perimental data p(3]. In the inset of Fig. 3 we plot d(u>) 
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using the parameters of the two cells of Fig.l, and both 
the high and low values of the effective surface tension (a 
and (To respectively) . Wc find a reasonable agreement be- 
tween the calculation and the measurements, when tak- 
ing D u ~ OAd. The similar magnitude of the bilayer- 
rigid shell separation from both static and dynamic ex- 
periments, shows the overall consistency of our confine- 
ment model. 

In the limit of high frequencies, the earlier result of 
Brochard et.al. 0, gave d(u>) 2 oc w~ 4 / 3 . Since qD M <C 1 
in the measured range, we find in this limit (|4|): uj(q) ~ 
nq^D^/Ar], leading to d{oj) 2 oc lu~ 7 / 5 . It is difficult to dis- 
tinguish between these two values using the newer data 
pof . Our calculation has the advantage of consistently 
describing both the static (spatial) and dynamic (tempo- 
ral) fluctuation data. Note that the case of a pure bilaycr 
with large effective surface tension er, but without the ef- 
fect of the rigid wall, is in complete disagreement with 
the data (dotted line, Fig. 3 inset). 

In Fig. 3 we show that the normal RBC, ATP depleted 
RBC and RBC ghost, are all well described by the same 
expression (|fj) (using the softer cell from Fig.l, i.e. the 
smaller 7 and a), differing by an effective temperature 
factor of up to ~ 3. This is similar to the amplitude 
enhancement factor of ~ 2.5 found in a previous study 
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The largest effect of the rigid shell is to increase the 
effective viscosity of the water near the bilayer, by con- 
straining its flow. Defining an effective viscosity by: 



Kq 3 /4T] eff , we get from (f|): ^ 
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the crossover wavevector go this function has its peak: 
Veff/v — 45 — 30, depending on the value of a. These 
values are in close agreement with the value rj e f f jr\ m 50 
found from the relaxation times of an electrodeformed 
RBC [H ||- In these experiments, the cytoplasm flows 
through the cytoskeleton mesh, setting up a flow field 



at the crossover wavevector qo. Thus, a rigid cytoskclc- 
tal wall, separated by a fixed distance from the bilaycr, 
accounts for the larger effective viscosity required to fit 
these dynamical experiments. 

While this model accounts for the wavevector depen- 
dence of the statics and the frequency dependence of the 
dynamics, the absolute amplitude of the fluctuations and 
the different values observed in active and ATP-depleted 
cells, must still be explained. One possibility is that 
ATP driven fluctuations completely determine the ampli- 
tude of the largest wavelength fluctuations |2(J through 



the process of s 
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spectrin-actin disconnections and recon- 
2p| . These ATP-driven conformational 
changes give rise to defects in the triangular spectrin net- 
work, resulting in nodes with more or less than 6 spec- 
trin bonds. The local curvature of the cytoskeleton may 
change at the site of a defect, from being locally flat (6 
bonds) to having a ~ 80nm deviation out of the plane of 
the flat cytoskeleton (5-fold node). The effect of this ran- 
dom buckling is to increase the mean bilayer-rigid shell 
separation by a factor of ~ 4. According to our model, 
this will increase the amplitude of the q — > modes by a 
factor of ~ 4 4 , as measured ||. 
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FIG. 1: The calculated (Eq. |2j) wavevector dependence of 
the bending modulus of the RBC (solid lines, taking 

a = sigmao) compared with the experimental data for the 
RBC P|(o,x). The crossover wavevector qo is indicated by the 
vertical dashed lines. Inset: A plot of (k 9 — K )~ 1 as a func- 
tion of the normalized wavevector (qd) 4 for small wavevectors. 
The linear slope in the limit of q — > is indicated by the dot- 
ted line. The deviation from linear behavior is well described 
by an effective surface tension a ~ k/£ 2 ~ 2.8, 7.7 x 10 -7 J/m 2 
for the two cells (solid line). Note that surface tension alone, 
without confining wall (7 = 0), does not describe the data 
(dash-dot line). 



5 




Normalized wavevector qd 



FIG. 2: A plot of the the measured effective modulus n/nq 
of RBC Q(o,x) and empty giant vesicles |15|(*) as a function 
of the normalized wavevector qd. The calculations are given 
by the solid and dashed line respectively. Inset: A plot of 
(K q — ft) -1 as a function of the normalized wavevector (qd) 4 
for small wavevectors. The linear slope in the limit of q — > 
for the RBC is indicated by the dotted line. The calculation 
for the vesicle is given by the solid line. 
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FIG. 3: Frequency dependence of the mean-square amplitude 
[ po[ d(io) 2 of the RBC (o), showing the reduction in the ampli- 
tude due to partial ATP depletion (x) and complete absence 
(RBC ghost)(*). The lines of the calculation (Eq. ||) differ 
by the effective temperature enhancement factor of ~ 3 (solid 
lines: do, dash-dot lines: a). Inset: A normalized Log-Log 
plot showing the powerlaw dependence (wo = 1Hz). The cal- 
culation is done using the parameters of the two cells of Fig.l 
(solid lines: ao, dash-dot lines: a). The dashed line shows 
the free bilayer behavior d(uj) 2 oc ui~ 5 ^ 3 . The case of a free 
bilayer with large effective surface tension a, but without the 
effect of the rigid wall, is in complete disagreement with the 
data (dotted line). 



